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ABSTRACT 

We present Spitzer images of the relatively sparse, low luminosity young clus- 
ter L988e, as well as complementary near-infrared (NIR) and submillimeter im- 
ages of the region. The cluster is asymmetric, with the western region of the 
cluster embedded within the molecular cloud, and the slightly less dense eastern 
region to the east of, and on the edge of, the molecular cloud. With these data, 
as well as with extant Ha data of stars primarily found in the eastern region 
of the cluster, and a molecular ^^CO gas emission map of the entire region, we 
investigate the distribution of forming young stars with respect to the cloud ma- 
terial, concentrating particularly on the differences and similarities between the 
exposed and embedded regions of the cluster. We also compare star formation in 
this region to that in denser, more luminous and more massive clusters already in- 
vestigated in our comprehensive multi-wavelength study of young clusters within 
1 kpc of the Sun. 
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1. INTRODUCTION 

Lynds 988e (L988e =IRAS 21023+5002) was identified as one of six IRAS point sources 
(L988a-f) observed toward the L98 8 dark cloud, which is located on the edge of the Cygnus 
0B7 molecular cloud association (jPobashi et al.l Il994l ) and which cont ains several bright, 
pre-main sequence objects. Bipolar molecular outfiows were discovered (jClarkill986l ) to be 
associated with the vicinity of L988e, as well as the vicinity of L988a and L988f, and there is 
also high velocity CO emission (blue-shifted lobe only) to the west of L988c. There are two 
LkHa objects and related nebulosity seen both at visible and infrared wavelengths, LkHa 
324-S E and LkHa 324 , the b rightest cluster members of a region of ~50 Ha emission-line 
stars (IHerbig fc Dahml |2006| ) on the edge of the heavy extinction associated with L988e's 
parent molecular cloud. Most of these emission-line stars lie to the east of the molecular 
cloud or on its eastern edge. LkHa 324 is a Herbig Ae/Be (HAeBe) star, close (in projection) 
to the IRAS coordinate, and LkH a 324-SE (misidentified as LkHa 3 24 in the the Herbig an d 
Bell catalog (IHerbig &: Belli Il988l )) is probably also an HAeBe star (IHerbig fc Dahml l2006l ). 
Both Classical T Tauri stars (CTTS) and Weak-line T Tauri stars (WTTS) are included 
in the emission-line stellar population. Employing an average extinction of Ay ~2.5 mag 



for rn ost of the population, and measured extinctions for a few objects, IHerbig fc Dahm 
( I2OO6I ) estimate median ages of 0.6 and 1.7 Myr for the emission-line stellar population 
from V and / photom e try: f or th ese estimat e s the y assume evolutionary isochrones from 
D'Antona fc Mazzitellil (119971 ). and lSiess et al.l (120001 ) respectively. In this paper, we examine 
properties of the entire cluster, including both the emission-line stellar population and an 
obscured population to the west within the molecular cloud which we have identified here, 
and refer to as the embedded population. 

At a distance of ~700 pc and a far infrared luminosity of ~ 200 L©, the cluster of 
young pre-main sequence stars associated with L988e is relatively sparse, and in the bottom 
quartile in far-infrared l uminosity of the n orthern young cluster sample surveyed in ^^CO 
and C^^O molecular gas (IRidge et al.l 120031 ) (see Table [ID. De spite present interes t in young 
cluster evolution as dis cussed, for exam ple, in a recent review (iLada fc Ladall2003l ). with the 
exception of papers by iHodappI (jl994l ) and IHerbig &: Dahm J2006h . the young cluster asso- 
ciated with L988e has been relatively unstudied at infrared wavelengths to date. Hereafter 
we refer to the young cluster as the L988e cluster, or just L988e. 

In this paper, we utilize mid-infrared (MIR) images of L988e obtained with the Infrared 
Array Camera (IRAC) and 24 /im images obtained with the Multi-Band Imaging Photome- 
ter (MIPS) on the Spitzer Space Telescop^ (see a composite three-color Spitzer image of 



^This work is based [in part] on observations made with the Spitzer Space Telescope, which is operated 
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the L988e region, also including L988f and c in Fig. [U also see Fig. [2] where the various 
fields of view are shown on a reproduction of the POSS V-band image of the region), as 
well as complementary near-infrared (NIR) ground-based images, to identify cluster mem- 
bers via their infrared excess emission from protostellar envelopes, and from circumstellar 
disks. The spatial distribution of identified protostars and disk objects in L988e is com- 
pared with the 850 /xm continuum maps of cold dust emission obtained with Submillimeter 
Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT), 
with dust extinction r naps generated us ing NIR data, and with a ^'^CO map of the associ- 
ated molecular clouds (IRidge et al.ll2003l ). These data were gathered as part of a coordinated 
multi-wavelength Young Stellar Cluster Survey, with the goal of providing the spatial dis- 
tribution of a complete census of YSOs - young stellar objects (all objects with an infrared 
excess, including stars with disks and protostars) in ~30 nearby e mbedded clusters chosen 
from a list of 73 young groups and clusters within 1 kpc of the Sun (jPorras et al.ll2003l ). Pro- 
gram clusters are representative of quite diverse star-forming regions which vary in number 
of cluster members, member number density, luminosity, and environment. This coordinated 
study, which includes NIR, mid-infrared (MIR), submillimeter and millimeter CO imaging 
of the regions, is designed to advance our understanding of clustered star formation. A com- 
prehensive paper examining YSO spatial distributions in this survey is in preparation: in a 
few instances, where unique data and/or characteristics distinguish the object, a separate 
paper is written. H ere, the availability of deep NIR images and a published catalog of Ha 
emission-line stars ( iHerbig &: Dahnj|2006l ) distinguishes this object from the others in the 
survey. 



Several authors have described efforts to define IRAC-only (lAllen et al.ll2004l ). or com- 
bined IRAC/MIPS ( iMuzeroUe et al.ll2004j ) protocols which differentiate between stars with 
dusty circumstellar accretion disks (class II), and protostars or objects dominated by dust 
emission from envelopes (class I). All class I and II objects exhibit infrared excess over 
photospheric emission typical of more evolved young cluster members (diskless or optically 
thin disk class III) as well as of foreground/background normal stars. Since substantial 
and/or variable interstellar extinction can cau se ambiguity in identification of these objects 
using IRAC-only methods JCutermuthl boosh . and the MIPS 24 /xm ir aages may not de- 
tect the faintest YSOs, a more robust classification method is required. iGutermuthI (120051 ) 
used seven-band photometry from J-band to IRAC 8.0 to systematically determine de- 
reddened SED slopes. Using this classification as a standard, methods requiring five-band 
detection at J (1.25 /im), H (1.65 /xm). Kg (2.15 /zm), [3.6] (3.6 /xm), [4.5] (4.5 /xm), or 
four-band detection at H ,Ks,[3.Q], [4.5] also have been demonstrated in the same publica- 



by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. 
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tion. Here we primarily employ a combination of the five-band and four-band methods, but 
we also use slopes deduced from IRAC photometry in combination with MIPS 24 fim pho- 
tometry, where available, to classify protostars and transition disk objects (objects with a 
pronour iced inner disk sa p) in L988e. Alternate methods are exploited for the young cluster 
IC 348 (ILada et al.ll2006l ): they utilize a simple power-law fit to the observed spectral energy 
distribution (SED) from the 4 IRAC bands to determine a slope ajRAc, and as well study 
the SED from visual to far-infrared wavelengths for all cluster members as a function of 
spectral class. 



2. OBSERVATIONS & DATA ANALYSIS 

2.1. Near-Infrared Data 

NIR images in the J, H, and Ks wavebands of L988e were first obtained by the authors 
on 2002 October 26 using the Simultaneous Quad Infrared Imaging Device (SQIID), which 
houses four 1024 x 1024 InSb arrays, of which only 512 x 512 are used, on the 2.1-m telescope 
at the Kitt Peak National Observatory, a facility of the National Optical Astronomical 
Observatory. The plate-scale of SQIID on the 2.1-m is 0".68 pixel"^, and the circular field of 
view is 5'. 8 in diameter. SQIID images, obtained in an overlapping 3x3 grid pattern, were 
used to construct a mosaic comparable in size to the 15' x 15' Spitzer IRAC field of view. 
Twelve dithered mosaics were obtained in each band with a total exposure time of 60 s per 
image at J, H, and Kg, leading to 720 s of exposure time. With seeing, the typical point 
spread function (PSF) was 1".9 FWHM at J-band and 1".77 in the other t wo bands, and the 



90% integrated completeness limits of the final mosaics were estimated (iGutermuth et al. 



20051 ) as Kg = 16.19, H = 16.20, and J = 17.71. However, the image quality suffered because 
all images exhibit position-dependent coma. PSF photometry was consequently impossible, 
and aperture photometry with large apertures was impractical because of the crowded field. 

Higher spatial resolution, deep, NIR observations in the J (1.25 /im), H (1.65 fim), 
and Kg (2.15 fim) wavebands were obtained on 2004 July 6 with the Palomar W ide-field IR 



Camera (WIRC) which houses a 2048 x 2048 HgCdTe array flWilson et al.ll2003l ). The plate- 
scale of WIRC is 0".249 pixel"^, and the field of view is 8'. 7 on a side. In addition, the PSF 
was l".l in the Kg band and 1".3 in the other two bands. The observing sequence consisted of 
a large (3'. 5 offset) five-point dither patte rn. The 90% integrated photometric completeness 



limits of the final mosaics were estimated (IGutermuth et al.ll2005l ) to be approximately Kg 
17.21, H = 17.72, and J = 18.22. 



Basic reduction of the SQIID NIR data was performed at Rochester using custom IDL 
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routines, developed by the authors, which include modules for linearization, fiat-field cre- 
ation and application, background frame creation and subtraction, distortion measurement 
and correction, and mosaicking. Point source detection and synthetic aperture photometry 
of all point sources wer e carried out using Phot Vis version 1.09, an IDL GUI-based photom- 



etry visualization tool (iGutermuth et al.ll2004l ). By visual inspection, those detections that 



were identified as structured nebulosity were considered non-stellar and rejected. Aperture 
photometry was performed using radii of 2". 04, 3". 40, and 5". 44 for the aperture, inner, and 
outer sky annular limits, respectively. 

WIRC data reduction was completed at Cornell using custom IDL routines and in- 
house image processing software. Processing steps included dark subtraction, linearization, 
sky subtraction, fiat-fielding, spatial fiux gradient correction (due to a dependency of the 
photometric zero point on position on the array; Jarrett personal comm.), and mosaicking. 
Aperture photometry using APPHOT was apphed to all point sources. Radii of the apertures 
and inner and outer limits of the sky annuli were 1".5, 2". 75, and 3". 5 respectively. PSF fitting 
was attempted, but the PSF varied sharply across the images, with resulting residuals larger 
than obtained for the aperture photometry discussed above. Aperture photometry with 
larger apertures was tested, but not applied because it failed in the crowded field. WIRC 
photometry was calibrated by minimizing residuals to corresponding 2MASS detections for 
selected stars with magnitudes between ~12-14 mag. Outliers were rejected from the fit 
(outliers are confused or unresolved multiples in the 2MASS dataset). The rms scatter for 



the J , H, and Kg waveband residuals are 0.06, 0.07, and 0.07 mag respectively. ICarpenter 



( 120011 ) showed that ~7% of YSOs toward Orion A are variable: variability may contribute to 
the observed scatter. Photometric magnitudes for stars which were brighter than the range of 
magnitudes corrected by the linearization modules were replaced with 2MASS photometry. 
Because of the superior spatial resolution, depth, image quality, yet similar field of view of the 
WIRC dataset, we use it for classification for all sources in L988e, and include photometry 
in Tables [21 [31 IH and [3 The SQIID dataset was used only for internal verification, and 
is not tabulated here. For fainter magnitudes (14-16) after rejecting outliers beyond 3(T, the 
rms scatter of the J, H, and Ks waveband residuals between the SQIID and WIRC datasets 
was 0.08, 0.11, and 0.11 mag respectively. We believe the WIRC photometry to be more 
accurate for the reasons cited above. 

Fig. [3] shows a three-co lor (WIRC J, iJ, Kg) image of the region with ^'^CO line emission 
contours (jllidge et al. 2003) overplotted. Visual inspection of the figure shows that the stellar 



density to the east of the molecular cloud is considerably higher than within the extent of 
the molecular cloud, as expected. 
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2.2. Mid-Infrared Data 

Mid-infrared images were obtained with Spitzer/IUKC at 3.6, 4.5, 5.8, and 8.0 /im and 
Spitzer /MIP'^ at 24 /xm as part of the Spitzer Young Stellar Cluster Survey. The IRAC 
images were obtained 2004 July 27, and each image was 5' on a side with 1".2 pixels. They 
were taken in 12 s High Dynamic Range (HDR) mode, whereby a 0.4 s integration is taken, 
immediately followed by a 10.4 s integration at the same pointing. Four dithered 3x4 
mosaics were obtained at each wavelength. The typical PSF is approximately 2". 2 at 3.6 
and 4.5 /im, 2". 3 at 5.8 fim, and 2". 5 at 8.0 fim. Ninety percent completeness limits of the 
final mosaics are approximately [3.6] = 14.68, [4.5] = 14.67, [5.8] = 13.72, and [8.0] = 12.73, 
and were evaluated using the same estimation technique as used for the NIR data. The 
Spitzer/MIPS 24 /im images were obtained on 2003 December 29, using scan mode at a 
medium scan rate. The total area imaged is about 20' x 1°. The typical point spread 
function is approximately 6".l, and the 90% completeness limit is [24] = 8.18. 

Basic reduction of the Spitzer/IRAC data was executed using the Spitzer Science Cen- 
ter's Basic Calibrated Data (BCD) pipelines and Spitzer/MIPS data using the Post-BCD 
product. The IRAC and MIPS Data handbook^ give complete descriptions of the calibration 
processes for point sources. The IRAC Data Handbook performs calibration photometry on 
a set of isolated standard stars using a 10 pixel (12") aperture, and 10 and 20 pixel radii of 
the sky annulus, but notes that considerably smaller apertures are required in crowded fields, 
and that small off-source sky annuli are needed when the point sources are embedded in ex- 
tended emission. In those cases, aperture corrections are required, and the Data Handbook 
3.0 provides aperture corrections for various standard aperture and annuli configurations for 
pipeline version S13. The same considerations apply to the MIPS data. 

Post-pipeline processing of IRAC data was performed at Rochester using custom IDL 
routines, developed by the authors, which include modules for bright source artifact correc- 
tion, cosmic ray removal and mosaicking. Point source extraction and aperture photometry 
were performed using PhotVis version 1.09. Radii of 2". 4, 2". 4, and 7". 2 were used for the 
aperture, inner, and outer sky annular limits, respectively, for the IRAC data and radii of 
7". 4, 14". 7, and 24". 5 were used for the aperture, inner and outer sky annular limits, respec- 
tively, for the MIPS data. For a Gaussian PSF, an aperture radius equal to the FWHM 
contains 93.7% of the total fiux. Increasing this aperture would negligibly gain signal from 
the source and increases total noise from the substantial MIR background emission. The 
two short-wave IRAC detector channels, 3.6 and 4.5 /im respectively, have approximately 



^The IRAC and MIPS Data Handbooks can be found at 
http: / / ssc. spitzer ■ caltech. edu / mips /dh/[ respectively. 



http: / / ssc. spitzer . caltech. edu /irac/dh/ 



and 
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Gaussian PSFs (other than the effects of diffraction), so that choosing an aperture radius 
approximately equal to the observed FWHM implies only a small correction because of the 
small aperture. Point source images for the two long-wave IRAC detector channels, 5.8 and 
8.0 nm respectively, do not have a Gaussian PSF: much of the point source intensity is scat- 
tered in bands along rows arid colu mns, as well as over the entire a rray and beyond (se e Data 
Handbook and lPipher et al.l (120041 )). but is well-characterized (see lReach et al.l (120051 )). The 
chosen aperture of 2". 4 only gives 63% of the source flux from a 12 " aperture in the 8.0 
/xm band. In this band, PAH emission nebulosity is very bright and often structured, neces- 
sitating use of relatively small apertures and sky annuli. IRAC photometry was calibrated 
using zero-magnitude flux densities of 280.9 ± 4.1, 179 7 ± 2.6, 115.0 ± 1 7, and 64.13 ± 0.94 
Jy in the 3.6, 4.5, 5.8, and 8.0 /im bands respectively (iReach et al.ll2005l ) in a 10 pixel (12") 
aperture, and aperture correction factors of 1.213, 1.234, 1.578, and 1.584 respectively as 
tabulated in the Data Handb ook. MIPS photometry was calibrated using a zero-point flux 
of 7.3 Jy, referenced to Vega JmIPS Data Handbook|[2005h . 



2.3. Submillimeter Data 

L988e was observed with the SCUBA submillimeter camera at 850 fim on the James 
Clerk Maxwell Telescop^ (JCMT) on Mauna Kea, Hawaii over the course of 5 nights from 
March to August 2003 in relatively poor weather conditions (r(225 GHz) = 0.25-0.5). The 
data were obtained in scan-mapping mode whereby a map is built up by sweeping the 
telescope across the source while chopping the secondary at a series of throws, 30", 44", and 
68". To minimize striping artifacts, each set of three sweeps was performed alte rnately in 



right ascension and declination. The SURF reduction package (jJenness et al.lll998l ) was used 
to flat-field the data, remove bad pixels, and make images. Although effective at removing the 
high sky and instrument background, this observing technique necessarily loses all features 
on scales larger than the largest chopper throw, 68". The calibration was performed by using 
the same observing technique to observe one or more of Mars, Uranus, and standard source 
CRL 2688. Fig. jH presents the SCUBA map in gray-scale, with ^^CO contours and YSOs 
(see below) overplotted. This map is presented here for morphological comparison to cluster 
member spatial distributions. 



■^The James Clerk Maxwell Telescope is operated by The Joint Astronomy Centre on behalf of the Particle 
Physics and Astronomy Research Council of the United Kingdom, the Netherlands Organisation for Scientific 
Research, and the National Research Council of Canada. 



- 8 - 



2.4. Extinction Maps 



We derive the Kks extinction maps using the {H — Kg) colors of all stars (11,313) 



detected in both bandpasses, using a method outlined in more detail by iGutermuth et al. 



(120051 ) ■ Briefly, the line of sight extinction to each point on a 2" grid throughout the region 
surrounding the L988e core was estimated using a variation of the nearest neighbor technique. 
At each point the mean and standard deviation of the {H - Kg) colors of the = 20 nearest 
stars to each grid point were calculated, using an iterative outlier rejection algorithm that 
would reject any star with an {H — Kg) color > 3 a from the mean, unt il the mean converged. 
The assumed NIR portion of the reddening law (IFlaherty et al.ll2007l ) given in Table E]was 
used to convert the {H — Kg) value at each point to A^^^, using A^s = 1.82((_f/' — Kg))obs - 
{H — Kg))int), where {H — Kg)int = 0.2 is assumed as an average intrinsic {H — Kg) color for 
stars in young clusters. Further, the H - Kg map was convolved with a 15" Gaussian kernel 
to match the SCUBA beara size. We compare the resultant map of Aks with the ^^CO 
map from iRidge et al.l (120031 ) in Fig. O As can be seen, the WIRC-derived A^^^ extinction 
map follows the general distribution of the molecular cloud as outlined by CO emission, but 
provides considerably greater spatial detail and extent. 

If stars within the molecular cloud are included in the above analysis, a biased estimation 
of Aj^g could result. In regions where cluster members dominate the A^=20 number counts, 
this estimation will be lower than the line of sight extinction through the entire cloud, 
and thsi deviation will be larger in regions of high extinction where stars beyond the cloud 
cannot be detected. To account for this bias, the {H — Kg) derived extinction map has 
been augmented with an extinction map derived by estimating A^g using SCUBA 850 /xm 
flux, and utilizing a conver sion factor of A/^-.s/Fs^n =0.272 mag pixel mJy~^ and an assumed 
dust temperature of 30 K (IGutermuth et al.ll2005l ). In areas where the SCUBA derived Kks 
values exceed the {H — Kg) - derived Aks values by 30%, the SCUBA derived values are 
used. Only one small region coincident with the L988e molecular core required augmentation 
with SCUBA-derived extinction values in the WIRC extinction map. 



3. IDENTIFYING YOUNG STELLAR OBJECTS 



3.1. NIR and MIR YSO identification Methods 



GutermuthI (120051 ) has shown that by de-redd ening the NIR/IR AC colors of young 



cluster stars, the IRAC-only classiflcation scher ae (lAUen et al.l l2004j) can be modifled to 
better isolate class II from class I objects. Briefly, GutermuthI ( 2005 ) reevaluated the original 
IRAC-only classiflcation scheme by considering the SED slopes from the NIR to the MIR. 
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Slopes (a = d log(AFx)/d log (A)) were constructed for both measured photometry and 
de- reddened photometry of point sources in the cluster field of view. When de-reddened 
photometry was utilized on six clusters, histograms of numbers of stars vs. 2-8 fim SED 
slope separate into 3 bins corresponding to field stars or class III (a < —2), class II (—2 < 
a < —0.5), and class I (a > —0.5), (see Fig. El). This figure is constructed using the 
reddening law discussed above (IFlaherty et al.l 120071 ). derived from IRAC and NIR cluster 
data: this reddening law proved m easurably different froni a previously published reddening 
law utilizing IRAC and NIR data (jindebetouw et al.ll2005l ) that is primarily associated with 
reddening from dust in the diffuse interstellar medium rather than from dust in a molecular 
cloud. 



In addition. iGutermuthI (120051 ) considered classification methods using J, H, i^rs,[3.6],[4.5] 
photometry as well as H, Kg, [3.6], [4.5] photometry. Where J-band d etections are availab le, 
objects are de-reddened to the (J - H)o vs. {H - Ks)o CTTS locus of iMeyer et al.l (119971 ) in 
the former method, or l ackin g a J-band det ection, to the {H - Ks)o vs. ([3.6] - [4.5] )o CTTS 
locus from IGutermuthI (120051 ) in the latter. ICutermuthl J2005h then applied the seven-band 
classifications to a de-reddened {K - [3.6] )o vs. ([3.6] - [4.5] )o color-color diagram, where the 
class I and II objects proved to be well separated. 

Since the 5.8 and 8.0 fim photometry is less sensitive, we utilize de-reddened (see TablelH 
for the assumed reddening law for young clusters) (J), H, Kg, [3.6], [4.5] photometry to 
classify most objects via their location on the {K — [3.6])o vs. ([3.6] — [4.5])o color-color 
diagram: these classifications are color-coded in Fig. 171 Where 24 fim data are available, 
objects classified as class I that have ([5.8] — [24]) < 4 were reclassified as class II. Additional 
class I and class II objects (hereafter class Is and class lis, similarly class Ills for class 
III objects) were identified utilizing SED slopes from IRAC and MIPS photometry, even 
when some objects were not detected in the H, Ks, or [3.6] bands. In addition, there are a 
handful of objects with NIR+IRAC colors of class Ills, but which exhibit a 24 /zm excess (e.g. 
([5.8] — [24]) > 1.5 and ([8] — [24]) > 1.5): these objects are classified as transition disk objects, 
YSOs characterized by a cool circumstellar disk containing a substantial inner gap. However, 
it must be noted that in the absence of a large 24 /xm excess, it is difficult to distinguish 
primordial transition disk objects (with ages less than a few Myr) from regenerative debris 
disk objects (with ages of tens of Myr to a few Gyr). Since the field of view covered by MIPS 
(as well as [3.6] and [5.8]) extends farther to the west than the (J), H, Kg, [3.6], [4.5] field, 
additional YSOs are identified outside that field. 
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3.2. Isolating Extragalactic Contamination 



Deep NIR and MIR photometric detections of YSOs toward molecular clouds are con- 
taminated by detection of extragalactic objects, including PAH (polycyclic aromatic hydro- 
carbon)-rich galaxies and active galactic nuclei (AGN), both of which have colors similar to 
YSOs. We account for extragalactic contamination in the following way. For objects which 
have detections in all four of the IRAC bands, we deredden each object on the assump- 
tion that it lies behind the molecular cloud, utilizing the extinction map discussed in §2.4. 
We then utilize ([4.5] - [5.8]) vs. ([5.8] - [8.0]) and ([3.6] - [5.8]) vs. ([4.5] - [8.0]) color- 
color diagrams to identify the PAH-rich galaxy contaminants and the [4.5] vs. ([4.5] — [8.0]) 
color-magnitude diagram to identi fy the AGN contaminant s. The prescriptions for identi- 
fication are described in detail by iGutermuth et al.l (120071 ) . For objects which only have 
(J)if/r[3.6][4.5] photometry, we are unable to directly use those prescriptions to identify 
extragalactic contaminants. Instead, inspection of the histogram of de-reddened [3.6] mag- 
nitudes for the objects with an infrared excess utilizing the ( J) if i^' [3. 6] [4.5] classification 
scheme (Fig. [S]) shows a bifurcation in the distribution at [3.6] = 13.5. Objects with [3.6] 
< 13.5 are considered as potential YSOs, while objects with [3.6] > 13.5 are considered as 
potential extragalactic contaminants and are rejected as potential YSOs. This scheme is in 
agreement with previous work stating that at magnitud es fainter than about 14, extragalac- 
tic objects dominate number counts (IFazio et al.ll2004l ). In order to quantify this further, 
we examine the differential num ber counts for sour ces detected in the Bootes region in the 
Flamingos Extragalactic Survey (jElston et al. 2006 ). At [3.6] = 13.5 (Fig. [8]), and a mean de- 
reddened {Ks — [3.6])o color of 0.79 (typical of YSO s), this bifurcation magnitude corresponds 
to Kg =14.29. From Fig. 8 in lElston et al.l (120061 ) . we estimate that there are approximately 
10 extragalactic sources within our field of view at Kg =14.29, and the number rapidly in- 
creases to fainter magnitudes. Since we are examining the spatial distribution as well as 
the disk fraction of the L988e population of YSOs, both in-cloud and off-cloud, we want to 
minimize the possibility of contamination by extragalactic objects. An examination of the 
Kg vs. {H — Kg) color-magnitude diagram (Fig. [H]) shows that even excluding Kg > 14.29, 
we are still capable of detecting relatively low mass YSOs. Faint cluster members (below the 
hydrogen burning limit at Ay=0) will also be excluded by these criteria: therefore brown 
dwarf cluster members will not be considered further. Spectroscopy is the best way to make 
firm brown dwarf identifications, and beyond the scope of this paper. 



3.3. Number Counts of Identified YSOs 



Before removing extragalactic contaminants, tliere are a total of 2293 objects in a 14.7' x 
12.5' field witli detections in all four bands (if, A'^, [3. 6], [4. 5]) meeting tlie criterion that the 
photometric uncertainties are less than 0.1 mag in all bandpasses. Following application of 
the extragalactic contamination criteria noted above, there remain a total of 678 objects in 
the NIR + IRAC dataset: these objects will be discussed in the following sections. Of these, 
three are class Is, 69 class lis, seven are transition disk objects, and 599 are class Ills or 
field stars with no obvious infrared excess. When SED slopes were computed for objects in 
the same field of view with only IRAC and/or MIPS detections, an additional five class Is, 
three class lis, and 36 class Ills were found. Further, there is one object detected in only 
the 4.5, 5.8, 8.0, and 24 /xm wavebands. From the 4.5 - 24 fim SED slope for this object, 
a=0.98, we characterize it as class I. As discussed earlier, a total of six transition disk objects 
are also classified based on their ([5.8] - [24]) or ([8] - [24]) colors. In Fig. [10] we present 
a three-color IRAC image ([3.6], [4.5], [8.0]) of the L988e region, with WIRC + IRAC class 
I and II classifications overplotted, and in Tables [21 [31 IH and \5\ we present a catalog of 
identified objects. 

The MIPS field of view extends far to the northwest of L988e, to include L988c, as noted 
above. Fortunately, IRAC data in two of the four bands (3.6 and 5.8 /im) were obtained 
to the northwest of L988e. From 3-24 fim SED slopes in combination with NIR data, a 
total of 19 class Is, 73 class lis, and 643 class Ills/field stars were found in the 14.7' x 19.4' 
field. A three-color IRAC + MIPS image ([3.6], [5.8], [24]) of the extended L988e region was 
presented in Fig. [TJ 



4. THE Ha EMISSION LINE PRE-MAIN SEQUENCE POPULATION 



Early studies of the emission-line stellar populations in molecular clo uds led to the 
classi cal T Tauri star (CTTS) and the weak T Tauri star (WTTS) definitions (iHerbig &: Bell 
19881 ) in terms of their Ha equivalent widths EW[IIa]. Although the WTTSs usually show 
weaker Ha emission, limited both in level and profile width (from chromospheric emission), 
CTTSs exhibit strong emission lines, with broad profiles, generated in the accretion process. 
By 1988, the CTTS class ification was ta k en to be EW[Ha] > 10 A. However, in a recent 
study of low mass stars, IWhite fc Basril (120031 ) suggest an empirical relationship for the 
critical EW[Ha] value as a function of spectral type for low mass stars. Using the presence 
of optical veiling to identify accretion, they propose the following improved scheme of CTTS 
classification for spectral classes later than KO: EW[Ha] > 3 A, for K0-K5; EW[Ha] > 
10 A, for K7-M2.5; EW[Ha] > 20 A, for M3-M5.5; and EW[Ha] > 40 A, for M6-M7.5. 
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In what follows, we first use the single equivalent width classification of CTTS/WTTS, and 
then examine whether this analysis would change appreciably with the lWhite fc Basril (120031 ) 
modification. 

In Fig. [5l cluster stars , including the HAeBe emission-line stars in this cluster identified 
by iHerbig fc Dahml (120061 ). are overplotted on an image of the extinction at Kg-hand (see 
§2.4). Each object in Fig. [5] is coded by YSO classification, deduced as outlined in §3 above: 
green squares (class II), red squares (class I), blue squares (transition disk objects). Filled 
squares denote YSOs that are Ha emission-line objects. Of course, with Ha emission line 
measurements, the degeneracy between field stars and pre-main sequence class HI cluster 
members is broken, so that class His with Ha emission can be identified (blue x's). The 
bright HAeBe stars LkHa 324-SE (V=15.36, Ay ~ 10) and LkHa 324 (V= 12.61, Ay = 
3.75) are identified by us as class I and class II, respectively. Even though the f^^-band error 
for LkHa 324 is slightly above our error threshold (error = 0.108, threshold = 0.1) we chose 
to include this object because of its prominence in the cluster. Most of the cluster associated 
with the emission-line stars lies in a region of low extinction to the east of the molecular 
cloud, or on the eastern edge of the cloud, where there is modest extinction, although there 
are bright, heavily extinguished members in the sample. 



In Fig.[TT], the Ha equivalent width EW[Ha], a proxy for accretion strength (IHerbig fc Dahm 



20061 ) ■ is plotted against the de-reddened SED slope fro m IRAC 3.6 to 8 . Atni photometry 
for the 41 emission-line stars common to the list from IHerbig fc Dahml (120061 ) and to our 
NIR -\- IRAC dataset with uncertainties < 0.1 mag in all four IRAC bands. A vertical line 
at EW[Ha] = 10 A divides the plot into the single-value estimate of W TTS (< 10 A) and 
CTTS (> 10 A). In their study of members of the young cluster IC 348, iLada et al.l (120061 ) 
show that diskless stars (Class HI) exhibit an IRAC determined SED slope of < -2.56, stars 
with optically thick disks exhibit slopes > -1.8, and they define sources with SED slopes 
between -1.8 and -2.56 to be 'anemic' disk objects (this class includes transition disk ob- 
jects, optically thin disks etc.): we have adopted these parameters in Fig. [TTJ Generally, 
Class lis (or Is) in L988e with o ptically thick disk s show the strongest Ha equivalent widths, 
as they do in IC 348 (see also iHartmannl (119981 ) who correlated {K — L) excess emission 
from optically thick disks with large EW[Ha]). In order to compare our results with earlier 
studies, we have indicated by color our Class I, II, and HI classifications. Although there is 
a general correlation of Class 11 /I with strong Ha equivalent widths in L988e, the figure also 
shows three stars with Ha equivalent widths > 10 A, but slopes indicating class HI accord- 
ing to the above definition, and three with Ha equivalent widths clearly < 10 A, but slopes 
indicating thick disks (Class II). There are five objects in the anemic disk region. Since six 
objects do not obey the correlation described above, we now determine whether use of the 
single- valued Ha equivalent width definition of CTTS/WTTS is the cause of the departure. 
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In the absence of spectroscopy for these 41 L988e cluster members, we made rough estimates 
of th e emission hne stars' spectral types from their measured {V-I) colors (IHerbig &: Dahm 
20061 ) for objects where good V and / photometry was available. In order to do so, we first 
confirm that spectral types determined from V and / photometry can be roughly calibrated 
for stars with disks in young clusters. 



The iLada et al.l (120061 ) study of cluster members in IC 348 included published spectro- 
scopic data, and they provide a table of median star-disk spectral energy distributions at 
wavelengths from l^-band to 24 /im, all with a reddening of Ay ~ 2.5, as a function of nine 
spectral type intervals. From their table, we confirm that the {V — I)o colors of both stars and 
objects with disks are sufficiently close for spectral typing into the nine interva ls. We also 



exam ined the JHKg colors for the same median star-disk population in IC 348 (ILada et al. 



20061 ) . and found NIR colors to be unsatisfactory for rough spectral typing of YSOs, because, 
as expected, those colors can vary widely for photospheres, stars with 'anemic disks' and 
stars with thick disks. 

Thirty-seven of the 41 emission line stars in L988e referred to above have good V and 
/ photometry. We first de-redden the V and I magnitudes using our determ ined values of 
Kks, and assuming Ky / Kks = 1 (values of this ratio can vary by up to 30% (lGlasslll999l )). 
Although iHerbig &: Dahml (120061 ) used a single value for Ay (2.5) with only a few exceptions. 



we found considerable variation in Kks for these stars. Thus {V — I)o colors were estimated 
for these 37 emission-line stars, and 62% are later than K7 as determined from those colors 
and the rough spectral type calibration. 



Applying the lWhite fc Basril (120031 ) criterion for CTTS/WTTS classification as a func- 
tion of rough spectral type, we found that 6 of the 37 objects yielded a different CTTS/WTTS 
classification than with the single- valued Ha equivalent width definition of CTTS/WTTS. 
These are indicated on Fig.[TT]with the appropriate symbol encircled, and color-coded (green 
circle implies CTTS is the correct classification rather than WTTS, and a blue circle implies 
WTTS rather than CTTS). The resultant correlation of Class III with WTTS and Class II 
with CTTS is superior to that obtained using the single-valued Ha equivalent width dis- 
criminator. Three transition disk objects are identified and are indicated by a blue square. 
One of these is an anemic disk and one is a thick disk object. The remaining transition disk 
object is found slightly below the -2.56 slope delineation, rather than above, so is marginally 
non-compliant with the correlation discussed above. There are two additional objects which 
fall in the anemic disk region. Overall, the correlation of CTTS/WTTS to Class Il/Class HI 
is excellent once the revision to the equivalent width CTTS/WTTS definition proposed by 
White fc Basril J2003h is applied. 



We now examine the statistics of the emission-line population. In addition to the 41 
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objects discussed above, there are four Ha objects which do not have detections in all four 
IRAC bands and six which do have detections, but with uncertainties in IRAC bands > 0.1 
mag. There are also two emission-line objects in common with our IRAC (3.6 and 5.8 fim) 
+ MIPS dataset in the vicinity of L988c. Of those objects with Ha emission, 28 are class lis, 
two are transition disk objects, 19 are class His, and two are class Is. Of the two in common 
with our IRAC (3.6 and 5.8 fim) + MIPS dataset, one is determined to be class I, and the 
other is determined to be class II from their 3-24 //m SED slopes and ([5.8] — [24]) colors. 

From these results, the disk fraction of just the emission-line population with NIR -|- 
IRAC counterpar ts is 67% ±20% including five which could be described as 'anemic disks' 
(ILada et al.ll2006l ) based on IRAC SED slopes or MIPS 24 /xm detections. This p ercentage 
is consistent with either of the two age estimates cited by iHerbig fc Dahml (I20061). namelv 
0.6 and 1.7 Myr, using the curve of disk fraction to age from [Hernandez et al.l (120071 ). 



5. THE COOL DUST DISTRIBUTION 

Young cluster environments are permeated by dust associated with the natal molecu- 
lar cloud, providing extinction as well as the cool dust emission detected at submillimeter 
wavelengths. All of the figures describing the molecular cloud will be used in the following 
discussion of the spatial distribution of the identified YSOs (Figs. [31111 andlSj). The dust and 
molecular gas are distributed non-uniformly, and in order to characterize that distribution 
for comparison with the YSO distribution, we map the extinction (Aks) in the neighborhood 
of L988e. The peak of the SCUBA 850 /xm emission (Fig. H]) coincides with the peak inten- 
sity of ^^CO emission, although a NE extension of the submillimeter emission does not have 
a parallel in the low spatial resolution molecular map: since the intensity of the emission 
varies on scales larger than the map size ( 5') outside the peak, there is no information on the 
larger scale dust emission due to the iterative way the map is built, assuming that the edges 
chop into zero emission. Thus the SCUBA-measured emission is not ideal for illustrating the 
extent of the molecular cloud. Despite adequately mapping the e xtent of the cloud, available 



low spatial resolution (50" beam size) CO emission line maps (IRidge et al.l 120031 ) are also 
not ideal for comparing the dust distribution with the YSO distribution. In addition to 
inadequate spatial resolution, there are issues associated with low signal to noise in diffuse 
regions (C^^O), and optical depth issues leading to saturation in dense cores (^'^CO), as well 
as CO depletion due to chemistry and/or freeze-out in dense cores, and/or varying excitation 
temperature. 

Although we could in principle utilize the IRAS 100 /xm HIRES map of the dust emission, 
it is more useful to use the higher spatial resolution [H — i^s)-derived extinction maps 
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described in §2 to compare the extended dust distribution with the YSO distribution. Note 
that the map of A^^^ (Fig- E]), at the same spatial resolution as the SCUBA map, shows the 
NE extension from the core indicating higher dust density, as seen by SCUBA. In addition, 
there are many small regions of enhanced extinction, and the extent of the molecular cloud 
defined by the Aks map exceeds that suggested by the CO maps. All cloud tracers show an 
abrupt eastern edge to the cloud, which bifurcates the central cluster into the relatively high 
surface density population embedded within the cloud core, and the somewhat lower surface 
density population outside the cloud. Th is edge was likely produced by the two HAeBe 
stars: as stated by iHerbig fc Dahiru (120061 ). LkHa 324 and 324SE 'have excavated a cavity 
on that edge and cleared out the area to the east, thereby exposing a small cluster of Ha 
emission-line stars.' 



6. SPATIAL DISTRIBUTION OF THE YOUNG STELLAR OBJECTS 

Classification of YSOs was executed as described in §3. We now examine the spatial 
distribution of the class Is and lis so identified, first by overplotting their positions on a 
variety of images: namely, the IR images, the CO image, the cold dust emission images, and 
primarily, the extinction map. 

Although the overall sie//ar density in the composite NIR image (Fig. 2) is substantially 
lower toward the molecular cloud, the peak surface density of YSOs is found toward the 
molecular cloud core, as identified by all of the molecular cloud tracers. LkHa324 and 
LkHQ!324SE are associated with that high density region, which we call the 'embedded' 
population of the cluster. 

To illustrate this more clearly, YSO surface density distributions are formed for all 
YSOs, utilizing the method developed for NGC 7129, GGD12-15 , and IRAS20050+2720 , 



three clusters previously studied in the Young Cluster Project (iGutermuth et al.l l2005l ). 
Because there exists a wide range in stellar surface densities over a young cluster, a method 
which utilizes adaptive smoothing lengths allows both the highest and lowest stellar density 
regions to be probed. At each sample position in a uniform grid the projected radial 
distance to the A^th nearest star, r]\i{i,j) is measured, and the local surface density <j{i,j) = 
N/ {irrj^^i, j)) is calculated. Here the grid size is 3", and N = 5. A value of > 5 is chosen 



to retain sufficient spatial resolution. On the other hand (IGutermuthl l2005l ) . A^ must be as 
sm all as feasible to a l low sm all, high density sub-groups to be properly represented. Also 
Gutermuth et al. J2004h . The grid size selection of 3" was chosen to fully resolve all 



see 



structure at the smallest scales (i.e. highest densities) in the A^ = 5 density distribution. 
The surface density of YSOs is compared with a SCUBA gray-scale image in Fig. [T2]of the 
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cool dust ridge associated with the edge of the molecular cloud. The surface density contours 
reveal an asymmetric cluster, with the highest surface density region (peaking at ~ 1706 
stars pc~^), including mostly class lis, concentrated at the peak of the dust emission. LkHa 
324 is located on the northern edge of this region and LkHa 324 SE lies on the eastern edge. 
A second, less compact surface density region associated with the emission-line population 
(peaking at ~ 747 stars pc~^) is found to the east of the dust ridge: these are also primarily 
class II objects (note that class I II objects are not incl uded in this figure, whereas in the 



emission-line selected population (IHerbig fc Dahmll2006l ) overp lotted on the map of A^s in 



Fig. O both class II and III stars are included). This population falls in the cavity where 
the molecular cloud has been excavated, presumably by the bright HAeBe stars LkHa 324 
and 324 SE. We call this the 'exposed' cluster population. 

The exposed population of the asymmetric cluster may be somewhat older than the 
embedded population, as suggested by the relatively lower YSO stellar densities encountered, 
as well as the necessity that enough time has passed to sculpt the cloud in the vicinity of the 
exposed cluster. We examine this further in the next section by estimating disk fractions for 
the exposed and embedded populations in a circular region surrounding the surface density 
peaks, using standard methods. 

In addition to these concentrated distributions of YSOs, in both the excavated cavity, 
and in the dense molecular knot, a distributed halo of YSOs is found mainly in the diffuse 
portions of the associated molecular cloud, and another to the east of the cloud. These 
haloes include predominantly class lis, but the presence of class Is in the diffuse molecular 
cloud is important to note. A much larger study focusing on the distributed population of 
YSOs well beyond the clusters, particularly the class Is, is underway. There we hope to 
determine whether the class Is in the distributed population formed in situ, or were ejected 
from the clusters. 



7. DISK FRACTIONS IN THE EXPOSED AND EMBEDDED CLUSTER 

REGIONS 

We show in this section that ~55±16% of the objects in a 1.5' diameter circle centered 
at [2lM™4.5*, 50°14'32".7] are class lis: this region corresponds to the peak of the exposed 
population of the cluster (here we have purposely not included the emission-line objects 
associated with substantial extinction at the edge of the L988e dark cloud, or the extended 
diffuse halo population, in contrast to the treatment in §4). This disk fraction, while esti- 
mated more rigorously than in §4, entails only a relatively small number of stars because 
L988e is sparsely populated. 
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In contrast, the disk fraction in a similar sized region centered on the embedded popula- 
tion's peak density (centered at [21^3'"52.97^ 50°14'43".5]) is ~74±15%. At face value, the 
higher mean disk fraction indicates the embedded population is younger than the exposed 
population, however, because of small number statistics, the large uncertainties overlap, so 
we cannot be certain of the relative ages of the regions. 

Disk fractions were estimated in the following manner, i^^-band magnitude histograms 
(KMH) were created for 1.5' diameter regions around the peaks of the exposed and embedded 
populations and for three 3' diameter regions off cloud that do not contain any known 
protostars, stars with disks, or Ha emission objects (see Fig. |2] for the locations of these 
regions and Fig. [13] for the histograms). To determine the magnitude at which field star 
contamination, begins to dominate the cluster KMH, we consider the faintest magnitude bin 
in the cluster KMH greater than Satackground, assuming Poisson counting statistics for each 
magnitude bin in the KMHs. This limit occurs at 14*^ magnitude. The average i^^-band 
extinction for the exposed population, the off-cloud regions, and the embedded population 
are Aks = 0.37, Aks = 0.28, and Aks = 1-3, respectively, estimated from J,H,Ks, colors as 
above (§3). To account for the higher extinction toward the embedded population, the off- 
cloud sources are artificially extinguished to the mean extinction of the embedded population 
before subtracting the off-cloud KMH from the embedded population KMH. By subtracting 
the off-cloud KMH from the exposed population KMH for the magnitude bins less than 
the 14*'^ magnitude limit, we estimate 22 cluster members in the exposed population, and 
19 cluster members in the embedded population: these values are used in the disk fraction 
determination above. 



8. COMPARISON WITH OTHER CLUSTERS 



We restrict comparison of L988e to six young clusters JCutermuthl boosh which have 



been studied in many of the same ways as L988e, although L988e is the only cluster with 
exposed stars which have been studied for Ha emission. Of those six, two filamentary 
clusters, GGD12-15 and IRAS 20050-1-2720, have the following properties: the class Is follow 
closely the asymmetric distribution defined by the dust ridges in their respective molecular 
clouds, while the class lis are more dispersed along the filament. The two clusters widely 
vary in far-IR luminosity and cloud mass, and range from moderate to relatively high spatial 
density (~ 1200 and 3000 pc'^ peak respectively). GGD 12-15 has 102 YSOs, of which 
25% are class I, while IRAS 20050+2720 has 121 YSO members, of which 21% are class 
I: class I number counts may be biased toward lower numbers due to the lack of MIPS 
coverage when these regions were analyzed. Two other clusters, NGC 7129 and AFGL 490, 
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are more diffuse than the filamentary clusters, with little high density clustering and little 
asymmetric structure. There are a total of 63 and 98 YSOs respectively in NGC 7129 and 
AFGL 490, of which 17% and 13% are class I respectively, with the same bias toward lower 
class I number counts due to lack of MIPS coverage. YSO spatial distribution in the L988e 
exposed cluster is similar to that in the latter two clusters, especially to that of NGC 7129: 
in particular, a large population of class lis clustered in a cavity in the molecular cloud 
near a Be star LkHa;234, resembles the exposed population discussed here. However, L988e 
does exhibit some asymmetry, and there is a substantial number of class Is in the L988e 
embedded cluster population. The distribution of class Is resembles those in GGD12-15 
and IRAS 20050+2720: class Is are found in both the high density dusty region near LkHa 
324 and 324-SE, and throughout the rest of the molecular cloud. However, in contrast to 
GGD12-15 and IRAS 20050+2720, only 10% of the 79 class I and II objects identified by 
IRAC and MIPS in all of L988e are class I. Of these four clusters, only IRAS 20050+2720 
has comparable FIR luminosity to L988e: the other three are ~ 6 - 30 times more luminous. 
NGC 7129 is estimated to be ~ 3 Myr in age, while the filamentary clusters are thought to be 
younger. The final two young clusters studied in this way, CepA and MonR2, have an order 
of magnitude higher FIR luminosity, but are otherwise different from each other. MonR2 
exhibits a filamentary distribution of 49 class Is (20% of the total) similar to that observed for 
GGD12-15 and IRAS 20050+2720, whereas CepA is sparsely populated by 8 class Is (only 9% 
of total YSOs) while the 77 class lis are anti-associated with CO density enhancements. A 
comprehensive spatial distribution study of the entire Young Cluster Survey is in preparation. 



9. SUMMARY 

1. S|?zfeer photometry of the asymmetric L988e cluster in combination with new NIR 
photometry reveals the presence of an embedded population associated with the molecular 
cloud (the cloud is defined by CO emission, cool dust emission, and Kks maps of extinction). 
This is in addition to the previously studied exposed region of the cluster, associated with 
the optically studied Ha emission-line stars in a cloud cavity on the eastern edge, sculpted 
presumably by the HAeBe stars LkHa 324 and LkHa 324 SE. 

2. The disk fraction (class II/(class II + class HI) of ONLY the emission-line objects, 
primarily located in the exposed region, is ~ 67%. This is consistent wit h either of the two 



age estimates deduced from V and I photometry (IHerbig fc Dahmll2006l ). 



3. The embedded population exhibits larger surface density than the exposed population 
of the cluster, and is fully enshrouded within L988e's natal molecular cloud. This, as well 
as the higher mean disk fraction calculated in the conventional way for all cluster members 
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in the field (74% vs. 55%) points to the probabihty of a somewhat younger age for the 
embedded population as compared with the exposed population, although the uncertainties 
of these disk fractions overlap. The exposed population of the cluster lies primarily in a 
region where the extinction A^s < 0.5 mag. On the other hand, the embedded population 
suffers a peak extinction of Kks = 3, and extinctions Aks > 0.5 throughout the molecular 
cloud. 

4. In addition to the asymmetric cluster, there exists a diffuse halo population of YSOs 
surrounding both the exposed and embedded regions of the cluster. Although the halo is 
dominated by class lis, class Is are also found in the diffuse halo. Whether these halo objects 
are formed in situ is yet to be determined. 

5. Comparin g these population s to six regions in the Young Cluster Program studied in 



a similar manner (iGutermuthl l2005l ) . we confirm the conclusion that the spatial distribution 
within the cluster is not strongly coupled to either the FIR luminosity or the numbers of 
cluster members. The younger clusters more faithfully follow the cloud density structure in 
which the cluster is born, and exhibit higher surface density. Therefore, L988e is older than 
the youngest clusters in the survey because the Eastern part of the cloud has had time to be 
excavated. It is becoming increasingly obvious that class Is are found in the haloes around 
clusters in regions where the molecular column density is relatively low. 

6. Spitzer/IKAC 3.6/im, 5.8/im and Spitzer/MlPS 24/xm photometry reveal another 
embedded cluster associated with the object L988c. 
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Table 1. L988e Cluster Characteristics Compiled from the Literature. 



L988o 


Cluster Distance (pc): 


700=^ 


Molecular Cloud Mass'' (Mq): 


197 


Velocity Dispersion'^ {km s~^): 


1.85(0.79) 


Far IR Luminosity'^ (^0): 


206.1 



'^ iClarklllisi) 

''Cl**0 core mass llRidge et al.ll2003h 

"^Value in parentheses is converted from 
Ql^O FWHM lin e width values for the core 
l lRidge et al.l2003h to equivalent Maxwellian val- 
ues. 
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Table 2. Objects with Mid Infrared Excesses (extragalactic contaminants removed): Class I 



1^-^-2000 






DEC2OOO 




J ^ 


H ^ 


K ^ 


[■^ fil ^ 


f4 51 ^ 


f5 81 ^ 


fs nl ^ 




[94.1 a 




CHIRAC ^ 


21:02:50.50±0 


23 


50 


18:52.43±0 


06 








10.50±0.003 




9.66±0.004 




5 


44±0.011 




-0.41S 


21:02:52.00±0 


32 


50 


13:43.96±0 


19 








10.59±0.004 




9.48±0.004 




5 


59±0.014 




-0.48S 


21:02:52.28±0 


05 


50 


12:29.47±0 


12 








8.55±0.002 




7.24±0.002 




2 


55±0.003 




0.020s 


21;02:54.17±0 


25 


50 


12:27.46±0 


10 








12.49±0.007 




10.01±0.005 




4 


36±0.008 




0.9608 


21:02:56. 60±0 


14 


50 


11:59.22±0 


04 








10.55±0.005 




9.65±0.007 




5 


40±0.014 




-0.37s 


21:02:58. 06±0 


19 


50 


15:25.63±0 


17 








9.91±0.002 




7.82±0.002 




3 


57±0.004 




0.080S 


21:03:00.19±0 


25 


50 


21:56.59±0 


13 








11.80±0.004 




9.56±0.004 




4 


03±0.006 




0.800S 


21:03:01. 16it0 


01 


50 


12:29.25±0 


03 








11.60±0.004 




10.19it0.005 




5 


67±0.018 




-0.04s 


21:03:02.98±0 


15 


50 


12:16.36±0 


18 








12.60±0.007 




11.59±0.014 




7 


62±0.076 




-0.47s 


21:03:03. 24±0 


12 


50 


13:12.66±0 


06 








9.79±0.004 




8.60±0.041 




2 


95±0.014 




0.470S 


21:03:03. 34±0 


43 


50 


09:31. 94±0 


15 








12.92±0.007 


12.26±0.006 


11.61±0.012 


10.48±0.011 








-0.05 


21:03:06. 85±0 


17 


50 


12:05.20±0 


37 








11.77±0.005 


11.29±0.004 


10.60±0.008 


9.54±0.005 


6 


72±0.022 




-0.25 


21:03:13. 90±0 


21 


50 


12:46.78±0 


06 








11.37±0.004 


10.57±0.004 


10.03±0.006 


9.20±0.005 


6 


16±0.015 




-0.40 


21:03:22.92±0 


22 


50 


11:31. 25±0 


23 








16.33±0.066 


15.19±0.052 


14.06±0.063 


12.96±0.067 


9 


60±0.151 




1.040 


21:03:35.66±0 


24 


50 


15:25.77±0 


33 




16.57±0.066 15.59±0.039 


13.92±0.030 


13.28±0.017 


12.92±0.096 


11.93±0.049 


6 


52±0.026 


1.09 


-0.64 


21:03:53. 44±0 


42 


50 


14:47.55±0 


39 










8.68±0.008 


7.00±0.018 


5.81±0.028 


1 


35±0.042 




0.98'' 


21:03:56. 37±0 


02 


50 


16:03.27±0 


08 








13.59±0.013 


11.25±0.005 


9.90±0.007 


8.62±0.007 


3 


97±0.024 




2.720 


21:03:58. 13±0 


34 


50 


14:39.77±0 


50 


7.71±0.015 


6.51±0.003 


5.32±0.011 


4.42±0.006 


3.64±0.005 


1.96±0.001 


1.65±0.001 






1.14 


0.580 


21:04:00.97±0 


48 


50 


15:56.77±0 


08 14.09±0.005 13.08±0.002 12.23±0.001 


10.25±0.002 


9.50±0.003 


8.79±0.009 


7.79±0.005 


4 


13±0.027 


0.68 


-0.03 



IHo'* EW(Hq:)' 



^Magnitudes are de-reddened when extinction value Axs quoted. 
''Ax determined from method outlined in §5. 

Slope a determined from IRAC photometry 

Running number from lHerbig &: Dahml l|2006l V 
<=EW[H«] from lHerbig fc Dahml l|2006l 'l. 
sSlope a determined from [3.6], [5.8] and [24] photometry 
^Slope a determined from [4.5] to [24] photometry 



Table 3. Objects with Mid Infrared Excesses (extragalactic contaminants removed): Class II 





RA2000 






DEC2000 






J 


a 




H =■ 


K ^ 


[3 fil » 

[O.OJ 




[5 81 » 

[O.OJ 


18 01 ^ 




1241 


a 




OilRAC ° 


IHa 


1—1 VV 1 


21 


02 


55 


21it0 


05 


50 


11 


55 


22i0 


22 














9.06±0.006 




7.68±0.004 




4 


34±0 


008 




-0 


66 


624 


29. 


21 


02 


59 


28it0 


03 


50 


10 


30 


24it0 


13 














11.61±0.004 




10.70±0.007 




6 


94±0 


061 




-0 


62^ 






21 


03 


06 


34it0 


24 


50 


10 


04 


04it0 


08 














11.06±0.004 


10.48±0.004 


9.99±0.005 


9.59±0.005 


7 


87±0 


047 




-1 


16 






21 


03 


08 


09±0 


15 


50 


10 


01 


71±0 


20 














14.31±0.020 


13.92±0.029 


13.57±0.062 


13.39±0.093 










-1 


78 






21 


03 


20 


22±0 


49 


50 


15 


58 


QlitO 


43 


9.94±0.008 


9.34±0.001 


9.22±0.001 


8.49±0.002 


8.14±0.003 


7.60±0.003 


6.64±0.002 


3 


55±0 


005 


2.15 


-0 


70 






21 


03 


21 


78it0 


27 


50 


09 


41 


73±0 


10 


15 


06±0 


Oil 


14.36±0.005 


14.05±0.007 


13.38±0.013 


13.18±0.029 


13.16±0.053 


13.17±0.070 








0.53 


-2 


64 






21 


03 


23 


72it0 


17 


50 


17 


21 


65i0 


56 


10 


24±0 


088 


9.60±0.004 


9.40±0.001 


9.09±0.003 


8.73±0.003 


8.29±0.005 


7.72±0.004 


5 


94±0 


021 


3.56 


-1 


25 






21 


03 


24 


17it0 


32 


50 


11 


01 


Olio 


34 














12.81±0.006 


12.11±0.005 


11.65±0.015 


11.77±0.061 










-1 


65 






21 


03 


26 


47±0 


37 


50 


14 


28 


25±0 


60 


15 


06±0 


012 


14 


22±0.005 


13.68±0.004 


12.79±0.008 


12.44±0.008 


12.03±0.020 


11.52±0.026 


7 


76±0 


043 


0.87 


-1 


37 






21 


03 


28 


67±0 


78 


50 


12 


34 


94±0 


47 


14 


26±0 


064 


13 


36±0.009 


12.70±0.005 


12.31±0.017 


12.16±0.015 


12.06±0.037 


12.13±0.088 








1.81 


-2 


64 






21 


no 


30 


'7/1 _l_n 
74±U 


70 


50 


12 


An 
49 


04ili0 


54 


12 


52±0 


046 


11 


80±0.004 


11.46±0.002 


11.00±0.005 


10.73±0.005 


10.22±0.010 


9.66±0.008 


7 


45±0 


051 


2.31 


-1 


zo 






21 


03 


36 


20it0 


48 


50 


15 


00 


D7±U 


33 


13 


24±0 


005 


12 


50±0.002 


12.12±0.001 


11.38±0.004 


10.98±0.004 


10.51±0.009 


9.92±0.010 


7 


83±0 


072 


1.07 


-1 


1 A 
14 






21 


03 


36 


24±0 


81 


50 


11 


50 


14i0 


53 


16 


57±0 


035 


15 


72±0.016 


15.15±0.015 


14.25±0.012 


13.89±0.015 


13.41±0.052 


12.56±0.064 








0.47 


-0 


90 






21 


03 


36 


69it0 


45 


50 


16 


49 


22i0 


13 


14 


58±0 


035 


13 


76±0.007 


13.23±0.005 


12.36±0.007 


12.12±0.007 


11.65±0.018 


11.12±0.022 


8 


30±0 


080 


1.40 


-1 


38 






21 


03 


37 


44it0 


56 


50 


13 


36 


68±0 


37 








13 


99±0.028 


13.69±0.015 


13.24±0.017 


13.12±0.015 


12.73±0.066 


12.56±0.081 








2.31 


-2 


01 






21 


03 


39 


47it0 


32 


50 


15 


52 


84it0 


27 


12 


40±0 


001 


11 


67±0.019 


11.30±0.022 


10.69±0.002 


10.52±0.003 


10.07±0.006 


9.48±0.008 


6 


78±0 


059 


0.14 


-1 


41 


631 


55. 


21 


03 


39 


68±0 


21 


50 


12 


36 


28±0 


55 








12 


40±0.028 


11.98±0.009 


11.50±0.008 


11.28±0.009 


10.78±0.024 


10.43±0.054 








3.26 


-1 


55 






21 


03 


40 


70±0 


11 


50 


13 


54 


05±0 


05 


11 


19±0 


001 


10 


54±0.001 


10.31±0.022 


9.73±0.003 


9.48±0.003 


9.09±0.004 


8.57±0.004 


6 


03±0 


026 


1.18 


-1 


49 






21 


Oo 


43 


co_i_r» 
08±U 


27 


50 


15 


CO 

00 


23±U 


00 


11 


01±0 


001 


10 


16±0.017 


9.60±0.019 


8.25±0.002 


7.82±0.002 


7.35±0.002 


6.97±0.002 


4 


25±0 


012 


0.77 


-1 


06 


632 


39. 


21 


03 


A A 

44 


23it0 


18 


50 


16 


05 


OOiO 


53 


14 


52±0 


030 


13 


92±0.008 


13.92±0.010 


13.17±0.039 


12.98±0.066 












1.46 










21 


03 


44 


46i0 


34 


50 


15 


51 


63i0 


33 


12 


28±0 


002 


11 


60±0.001 


11.31±0.001 


10.66±0.009 


10.23±0.007 


9.81±0.010 


9.22±0.009 








0.99 


-1 


19 


633 


156 


21 


03 


46 


39±0 


24 


50 


13 


44 


79±0 


05 


15 


61±0 


003 


14 


77±0.002 


14.22±0.002 


13.50±0.015 


13.15±0.018 


12.81±0.050 


12.35±0.084 


7 


85±0 


155 


0.09 


-1 


52 






21 


03 


46 


58±0 


27 


50 


08 


41 


97±0 


82 


12 


51±0 


001 


11 


85±0.001 


11.60±0.001 


11.32±0.004 


11.15±0.004 


10.87±0.012 


10.25±0.036 


7 


57±0 


115 


0.15 


-1 


62 


635 


4.2 


21 


03 


46 


71±0 


43 


50 


15 


26 


33±0 


24 








13 


22±0.009 


12.33±0.004 


11.01±0.005 


10.43±0.004 


9.72±0.014 


8.98±0.035 


6 


47±0 


186 


1.87 


-0 


49 






21 


03 


47 


19±0 


20 


50 


11 


27 


63±0 


29 


13 


25±0 


001 


12 


62±0.001 


12.44±0.001 


12.31±0.006 


11.94±0.005 


11.61±0.015 


10.73±0.015 


8 


02±0 


116 


0.19 


-1 


06 


637 


19. 


21 


03 


48 


42±0 


37 


50 


14 


18 


15±0 


13 


11 


08±0 


oil 


10 


29±0.001 


9.84±0.001 


9.11±0.003 


8.68±0.003 


8.15±0.004 


7.44±0.003 


3 


32±0 


006 


1.97 


-0 


91 






21 


03 


48 


89±0 


12 


50 


10 


04 


72±0 


54 


13 


33±0 


001 


12 


64±0.001 


12.33±0.001 


11.96±0.005 


11.78±0.005 


11.65±0.019 


10.94±0.055 


6 


20±0 


043 


0.27 


-1 


71 


638 


28. 


21 


03 


49 


03±0 


14 


50 


16 


02 


59±0 


26 


12 


13±0 


002 


11 


25±0.001 


10.62±0.001 


9.55±0.002 


9.23±0.003 


8.97±0.004 


8.56±0.016 


6 


26±0 


064 


0.95 


-1 


72 


639 


6.8 


21 


03 


50 


14±0 


01 


50 


10 


23 


79±0 


06 


13 


26±0 


001 


12 


45±0.001 


11.95±0.001 


11.26±0.004 


10.85±0.003 


10.48±0.007 


9.85±0.010 


8 


10±0 


286 


0.22 


-1 


24 


640 


86. 


21 


03 


50 


80±0 


08 


50 


15 


33 


02±0 


19 


14 


12±0 


077 


13 


35±0.011 


12.92±0.006 


12.68±0.047 


12.50±0.040 












1.94 










21 


03 


51 


50±0 


11 


50 


14 


32 


72±0 


31 


10 


81±0 


103 


10 


18±0.006 


9.62±0.001 


8.68±0.005 


8.37±0.004 


7.99±0.016 


7.63±0.041 


3 


95±0 


021 


3.29 


-1 


62 






21 


03 


52 


17±0 


02 


50 


14 


43 


89±0 


34 


12 


40±0 


005 


11 


30±0.001 


10.30±0.001 


9.12±0.007 


8.58±0.006 


8.03±0.026 


7.54±0.059 








1.33 


-1 


00 






21 


03 


52 


89±0 


05 


50 


14 


38 


07±0 


09 


14 


42±0 


004 


13 


64±0.003 


13.18±0.004 


12.51±0.070 


12.06±0.078 












0.25 










21 


03 


53 


32±0 


31 


50 


17 


49 


99±0 


32 


16 


96±0 


030 


16 


20±0.018 


15.80±0.020 


14.66±0.034 


14.52±0.028 


13.88±0.078 


13.32±0.090 








0.44 


-1 


21 






21 


03 


53 


56±0 


31 


50 


16 


00 


64±0 


07 


14 


05±0 


003 


13 


24±0.002 


12.75±0.001 


12.07±0.006 


11.69±0.008 


11.51±0.045 


10.80±0.106 








0.47 


-1 


44 







to 



Table 3 — Continued 





RA2000 






DEC2000 




T a 
J 


IT a 
n 




[3.bJ 


[4.5J 


[5.SJ 


[S.UJ 








OtIRAC 


illCK 


tl W \ n.Oi 


21 


03 


53 


87±0 


01 


50 


14 


44 


48±0 


01 


11.50±0.006 


10.77±0.002 


10.41±0.001 


9.86±0.034 


9.49±0.029 










1 


44 








21 


03 


53 


91±0 


00 


50 


15 


14 


39±0 


23 


14.12±0.005 


13.32±0.004 


12.84±0.004 


11.63±0.062 


11.10±0.051 













28 




643 


43. 


21 


03 


54 


10±0 


05 


50 


15 


29 


26±0 


11 


12.24±0.002 


11.53±0.001 


11.21±0.001 


10.49±0.009 


10.00±0.007 


9.65±0.074 


8.70±0.158 






1 


06 


-0.84 






21 


03 


54 


21±0 


32 


50 


15 


10 


15±0 


13 


10.02±10.00 


9.81±0.015 


9.57±0.090 


9.51±0.017 


9.27±0.018 


8.50±0.109 











00 






11. 


21 


03 


54 


36±0 


16 


50 


15 


00 


07±0 


01 


12.57±0.002 


11.94±0.001 


11.74±0.001 


11.00±0.094 


10.64±0.056 













45 




644 


15. 


21 


03 


54 


43±0 


37 


50 


13 


50 


73±0 


28 


11.90±0.003 


11.26±0.001 


11.04it0.001 


10.47±0.004 


10.30±0.006 


10.10±0.019 


9.79it0.032 






1 


47 


-2.06 






21 


03 


54 


63±0 


11 


50 


14 


36 


87±0 


25 


8.79±0.025 


7.89±0.021 


7.23±0.037 


5.88±0.002 


5.25±0.002 


4.51±0.001 


3.56±0.003 





95±0.006 


1 


49 


-0.16 


645 




21 


03 


55 


05±0 


30 


50 


14 


58 


53±0 


12 


12.84±0.007 


11.97±0.003 


11.35±0.002 


10.09±0.046 


9.60±0.023 










1 


24 








21 


03 


55 


63±0 


04 


50 


14 


14 


28±0 


09 


15.21±0.024 


14.26±0.008 


13.51±0.006 


12.18±0.044 


11.75±0.071 


11.45±0.165 


10.14±0.142 









95 


-0.57 






21 


03 


56 


11±0 


01 


50 


10 


02 


70±0 


14 


13.78±0.002 


12.91±0.001 


12.31±0.001 


11.50±0.005 


11.24±0.004 


10.88±0.014 


10.69±0.047 









00 


-1.89 


647 


19. 


21 


03 


56 


33±0 


14 


50 


14 


51 


39±0 


20 


12.00±0.001 


11.21±0.001 


10.73±0.001 


9.53±0.007 


9.15±0.012 


8.58±0.038 


8.79±0.285 









24 




648 


31. 


21 


03 


56 


42±0 


27 


50 


10 


03 


46±0 


32 


12.88±0.001 


12.22±0.001 


11.98±0.001 


11.77±0.006 


11.49±0.005 













06 








21 


03 


56 


95it0 


57 


50 


13 


51 


61±0 


28 




13.72±0.045 


13.47it0.020 


12.90it0.046 


12.82±0.087 










2 


75 








21 


03 


56 


96±0 


08 


50 


15 


17 


71±0 


13 


13.40±0.001 


12.67±0.001 


12.30±0.001 


11.45±0.017 


10.84±0.015 


10.72±0.166 


10.12±0.562 









35 




649 


233 


21 


03 


58 


90±0 


13 


50 


14 


50 


88±0 


34 


11.33±0.001 


10.55±0.024 


10.12±0.001 


9.06±0.011 


8.54±0.029 


8.25±0.072 


7.82±0.157 









33 


-1.47 


652 


53. 


21 


03 


59 


24±0 


37 


50 


14 


57 


66±0 


15 


11.16±0.001 


10.48±0.001 


10.21±0.001 


9.47±0.012 


9.13±0.017 


8.60±0.134 


7.91±0.132 









75 


-1.02 


654 


53. 


21 


04 


00 


18±0 


36 


50 


15 


24 


36±0 


09 


12.47±0.001 


11.70±0.017 


11.28±0.022 


10.45±0.003 


10.02±0.003 


9.54±0.028 


8.82±0.040 









11 


-0.96 


655 


98. 


21 


04 


00 


94±0 


17 


50 


17 


08 


63it0 


36 


13.77±0.019 


13.06it0.005 


12.72±0.004 


12.44±0.011 


12.20±0.013 


12.16±0.044 


12.50±0.125 






1 


39 


-2.93 






21 


04 


00 


95±0 


20 


50 


17 


21 


39±0 


17 




12.36±0.011 


12.08±0.007 


11.66±0.039 


11.55±0.042 


10.62±0.053 








2 


40 








21 


04 


01 


02±0 


02 


50 


17 


27 


05±0 


08 


10.94±0.002 


9.95±0.001 


9.16±0.021 


8.12±0.003 


7.64±0.002 


7.18±0.002 


6.64±0.002 


4 


53±0.008 


1 


58 


-1.15 






21 


04 


01 


34±0 


32 


50 


14 


41 


85±0 


30 


12.63±0.001 


11.74±0.035 


11.01±0.031 


9.81±0.004 


9.42±0.005 


8.80±0.032 


8.28±0.041 









00 


-1.03 


658 


4.1 


21 


04 


01 


50±0 


08 


50 


15 


05 


13±0 


25 


13.52±0.002 


12.75±0.001 


12.30±0.001 


11.51±0.010 


11.08±0.020 


10.74it0.052 


9.86±0.076 









23 


-0.99 


659 


183 


21 


04 


01 


77±0 


11 


50 


14 


23 


22±0 


38 


13.44±0.001 


12.76±0.001 


12.47±0.001 


12.24±0.029 


12.05±0.036 


11.76±0.267 











27 








21 


04 


01 


87±0 


53 


50 


14 


43 


11±0 


27 


13.21±0.001 


12.40±0.001 


11.89±0.001 


11.25±0.012 


10.82±0.014 


10.29±0.077 


9.28±0.051 









26 


-0.58 


660 


121 


21 


04 


02 


03±0 


24 


50 


14 


10 


lOitO 


27 


14.92±0.003 


14.20±0.002 


13.86±0.003 


13.45±0.041 


13.17±0.055 













13 








21 


04 


02 


77±0 


19 


50 


14 


20 


31±0 


05 


11.08±0.013 


10.35±0.015 


9.97±0.019 


8.91±0.003 


8.36±0.002 


7.63±0.005 


6.66±0.010 


3 


27±0.027 





12 


-0.23 


663 


12. 


21 


04 


02 


80±0 


18 


50 


15 


29 


56±0 


24 


12.17±0.001 


11.55±0.001 


11.38±0.001 


11.02±0.005 


10.54±0.004 


10.00±0.012 


9.26±0.020 


5 


48±0.127 





65 


-0.81 


664 


9.6 


21 


04 


03 


36±0 


17 


50 


14 


44 


47±0 


12 


13.24±0.001 


12.34±0.021 


11.69±0.001 


10.77±0.004 


10.35±0.004 


9.83±0.024 


9.03±0.016 









03 


-0.84 


666 


46. 


21 


04 


05 


50±0 


16 


50 


14 


36 


15±0 


24 


15.54±0.003 


14.81±0.003 


14.35±0.003 


13.63±0.028 


13.16±0.025 


12.51±0.085 


12.00±0.046 









00 


-0.93 






21 


04 


05 


72±0 


38 


50 


14 


47 


85±0 


19 


13.30±0.002 


12.65±0.001 


12.44±0.001 


11.85±0.011 


11.59±0.007 


11.19±0.030 


10.56±0.056 


6 


86±0.162 





28 


-1.34 






21 


04 


06 


11±0 


09 


50 


16 


39 


01±0 


34 


12.62±0.001 


11.85±0.001 


11.42±0.001 


10.23±0.004 


9.83±0.003 


9.46±0.004 


8.62±0.006 


5 


45±0.017 





53 


-1.01 


670 


28. 


21 


04 


06 


82±0 


09 


50 


13 


17 


63±0 


34 


11.33±0.001 


10.65±0.015 


10.37±0.021 


10.11±0.003 


9.55±0.003 


9.11±0.004 


8.25±0.004 


5 


12±0.031 





73 


-0.74 


671 


161 


21 


04 


06 


86±0 


09 


50 


14 


41 


69±0 


16 


13.72±0.001 


12.93±0.001 


12.48±0.002 


11.78±0.009 


11.43±0.007 


11.04±0.022 


10.51±0.060 


7 


44±0.173 





04 


-1.38 


672 


101 


21 


04 


07 


05±0 


20 


50 


18 


52 


32±0 


12 


12.01±0.004 


11.41±0.001 


11.29±0.001 


10.76±0.004 


10.65±0.004 


10.15±0.010 


9.24±0.006 


5 


57±0.013 


1 


32 


-1.04 






21 


04 


07 


09±0 


05 


50 


15 


19 


01±0 


17 


12.96±0.002 


12.35±0.002 


12.20±0.002 


11.80±0.009 


11.55±0.012 


11.18±0.014 


10.49±0.026 


6 


92±0.059 





41 


-1.33 


673 


17. 



Table 3 — Continued 



RA2000 



DEC2000 



H 



K 



[3.6] 



[4.5] 



[5.8] 



[8.0] 



[24] 



Aks'' a/flAc" IHa<i EW(Ha) ' 



21:04:11. 23±0.03 50:10: 
21:04:19.50±0.24 50:15: 
21:04:20.57±0.07 50:16: 
21:04:28.81±0.39 50:18: 



48.57±0.07 12.03±0.001 11.33±0.017 11.03±0.021 10.21±0.002 9.94±0.003 9.56±0.004 8.94±0.004 6.97±0.031 0.27 -1.38 

57.86±0.22 11.88±0.001 11.00±0.001 10.40±0.019 9.39±0.002 9.14±0.003 8.73±0.003 8.04±0.002 5.21±0.010 0.05 -1.27 

08.26±0.13 13.83±0.025 12.96±0.006 12.35±0.004 11.14±0.005 10.73±0.007 10.08±0.009 9.20±0.008 6.37±0.027 1.41 -0.57 

24.20±0.07 13.48±0.002 12.70±0.001 12.25±0.001 11.61±0.004 11.38±0.004 10.98±0.008 10.14±0.008 7.44±0.036 0.20 -1.15 



676 
681 



26. 
35. 



''Magnitudes are de-reddened when extinction value Aks quoted, 
determined from method outlined in §5. 

Slope a determined from IRAC photometry 

Running number from iHerbig fc Dahml l l2006h . 
''EWiHo] from lHerbig fc Dahml l|2006l '). 
'slope a determined from [3.6], [5.8] and [24] photometry 
^Slope a determined from K-band to longest detected wavelength 



Table 4. Objects with Mid Infrared Excesses (extragalactic contaminants removed): Transition Disks 



RA2000 DEC2000 J " H ^ K ^ [3.6] ^ [4.5] ^ [5.8] ^ [8.0] ^ [24] Aa'. a/flAC " IHa ^ EW(Ha) < 

21;03:27.22±0.76 50:13;00.70±0.41 11.63±0.001 10.97±0.001 10.72±0.024 10.65±0.003 10.52±0.003 10.35±0.008 9.80±0.006 6.60±0.022 0.88 -1.88 

21:03:52.07±0.18 50:14;18.65±0.36 11.19±0.001 10.59±0.017 10.55±0.021 10.62±0.005 10.71±0.009 10.62±0.030 10.49±0.056 4.83±0.027 0.59 -2.67 641 17 

21;03:57.80±0.27 50:12:31. 18±0.21 14.00±0.003 13.37±0.001 13.17±0.002 12.94±0.007 12.83±0.008 12.82±0.046 12.53±0.083 7.92±0.107 0.38 -2.40 650 37 

21:04:02. 41±0.44 50:07:39. 27±0.58 6.58±0.011 5.95±0.015 5.75±0.007 5.64±0.002 5.86±0.002 5.52±0.002 5.45±0.005 3.84±0.018 0.85 -2.52 

21:04:08. 69±0.03 50:14:25. 80±0.20 12.79±0.001 12.19±0.001 12.07±0.001 11.79±0.014 11.64±0.010 11.38±0.073 10.56±0.048 8.17±0.225 0.36 -1.43 674 12 

21:04:12. 64±0.01 50:16:00.98±0.24 14.56±0.007 13.59±0.002 12.80±0.002 11.46±0.006 11.50±0.008 10.54±0.008 10.19±0.014 7.83±0.068 0.48 -1.17 

21:04:23.84±0.46 50:13:05. 78±0.40 14.06±0.002 13.35±0.001 13.03±0.002 12.53±0.007 12.49±0.007 12.30±0.021 11.98±0.048 9.29±0.129 0.03 -2.19 

''Magnitudes are de-reddened when extinction value A^s quoted. 
^Ajf determined from method outhned in §5. 

Slope a determined from IRAC photometry 
•^Running number from lHerbig &: Dahml ll2006ll . 
°EW[Ha] from lHerbig fc Dahml ||2006|V 



Table 5. Objects without Mid Infrared Excesses (extragalactic contaminants removed): Class III 





RA2000 






DEC2000 






J 


a 




}{ 


a 






a 




f3 61 


a 




f4 51 


a 




f5 81 


a 




fS 01 

[O.UJ 


a [24.1 a 




aiRAC " 


IHa 


EW(Ha) ° 


21 


03 


29 


19ib0 


11 


50 


19 


05 


84ib0 


26 


14 


48±0 


003 


13 


88±0 


003 


13 


92±0 


005 


13 


79±0 


013 


13 


85±0 


018 


13 


94±0 


091 











29 






628 


1.3 




21 


03 


32 


1 / ±U 


68 


50 


10 


55 


25i0 


63 


14 


52±0 


002 


13 


92±0 


001 


13 


66±0 


002 


13 


19±0 


008 


13 


15±0 


010 


13 


10±0 


042 


13 


05±0 


073 





00 


-2 


68 


629 


3.9 




21 


03 


33 


94±0 


26 


50 


20 


32 


81±0 


00 


14 


25±0 


002 


13 


65±0 


002 


13 


57±0 


003 


13 


35±0 


010 


13 


29±0 


Oil 


13 


12±0 


042 


13 


16±0 


078 





12 


-2 


60 


630 


3.5 




21 


03 


46 


60i0 


27 


50 


10 


43 


50i0 


25 


12 


36±0 


001 


11 


69±0 


001 


11 


43±0 


001 


11 


19±0 


005 


11 


19±0 


004 


11 


12±0 


014 


11 


14±0 


037 





50 


-2 


( / 


636 


11. 




21 


03 


52 


07i0 


18 


50 


14 


18 


65i0 


36 


11 


19±0 


001 


10 


59±0 


017 


10 


55±0 


021 


10 


62±0 


005 


10 


71±0 


009 


10 


62±0 


030 


10 


49±0 


056 4.83±0.027 





59 


-2 


67 


641 


17. 




21 


03 


57 


80i0 


27 


50 


12 


31 


18i0 


21 


14 


00±0 


003 


13 


37±0 


001 


13 


17±0 


002 


12 


94±0 


007 


12 


83±0 


008 


12 


82±0 


046 


12 


53±0 


083 7.92±0.107 





38 


-2 


40 


650 


37. 




21 


04 


00 


27±0 


15 


50 


14 


20 


09±0 


15 


12 


85±0 


001 


12 


23±0 


001 


12 


05±0 


001 


11 


84±0 


010 


11 


82±0 


035 


11 


65±0 


138 











29 






656 


9.5 




21 


04 


02 


57±0 


06 


50 


14 


08 


70±0 


40 


12 


45±0 


001 


11 


85±0 


001 


11 


73±0 


001 


11 


63±0 


015 


11 


60±0 


026 


11 


36±0 


063 


11 


22±0 


156 





39 


-2 


33 


662 


4.5 




21 


04 


03 


69±0 


14 


50 


13 


23 


04±0 


01 


12 


02±0 


001 


11 


38±0 


019 


11 


17±0 


024 


11 


02±0 


004 


10 


92±0 


004 


10 


85±0 


012 


10 


82±0 


021 





25 


-2 


61 


665 


19. 




21 


04 


04 


16±0 


09 


50 


15 


01 


24±0 


20 


11 


37±0 


001 


10 


77±0 


017 


10 


67±0 


021 


10 


55±0 


003 


10 


57±0 


004 


10 


46±0 


016 


10 


53±0 


048 





50 


-2 


80 


667 


1.8 




21 


04 


04 


50±0 


03 


50 


14 


01 


54±0 


09 


12 


97±0 


001 


12 


32±0 


001 


12 


08±0 


001 


11 


73±0 


005 


11 


71±0 


006 


11 


62±0 


036 


11 


49±0 


100 





21 


-2 


56 


668 


8.6 




21 


04 


05 


58±0 


01 


50 


15 


58 


35±0 


10 


12 


04±0 


001 


11 


44±0 


060 


11 


29±0 


024 


11 


07±0 


005 


11 


15±0 


006 


11 


01±0 


013 


11 


05±0 


034 





39 


-2 


78 


669 


8.6 




21 


04 


08 


80±0 


26 


50 


14 


33 


10±0 


05 


12 


76±0 


001 


12 


14±0 


001 


11 


96±0 


001 


11 


69±0 


008 


11 


71±0 


009 


11 


62±0 


039 


11 


65±0 


050 





22 


-2 


77 


675 


2.7 




21 


04 


14 


19±0 


05 


50 


09 


52 


96±0 


10 


13 


28±0 


001 


12 


68±0 


001 


12 


60±0 


001 


12 


37±0 


006 


12 


38±0 


006 


12 


31±0 


022 


12 


10±0 


037 





07 


-2 


53 


677 


5.1 




21 


04 


16 


19±0 


05 


50 


13 


34 


94±0 


04 


13 


48±0 


001 


12 


88±0 


001 


12 


76±0 


001 


12 


59±0 


010 


12 


66±0 


010 


12 


56±0 


036 


12 


61±0 


046 





27 


-2 


83 


678 


4.3 


to 

CO 


21 


04 


18 


92±0 


09 


50 


20 


28 


47±0 


39 


12 


33±0 


002 


11 


73±0 


001 


11 


69±0 


001 


11 


37±0 


004 


11 


61±0 


006 


11 


28±0 


017 


11 


12±0 


023 





74 


-2 


44 


679 


35. 


21 


04 


19 


18±0 


05 


50 


13 


18 


06±0 


04 


12 


89±0 


001 


12 


27±0 


001 


12 


03±0 


001 


11 


75±0 


004 


11 


65±0 


004 


11 


60±0 


014 


11 


42±0 


032 





00 


-2 


48 


680 


2.1 




21 


04 


22 


07±0 


08 


50 


15 


44 


40±0 


20 


11 


68±0 


001 


11 


08±0 


017 


10 


95±0 


022 


10 


85±0 


004 


11 


00±0 


004 


10 


88±0 


Oil 


10 


75±0 


013 





23 


-2 


69 


682 


2.1 




21 


04 


23 


23±0 


20 


50 


13 


00 


74±0 


34 


12 


63±0 


001 


12 


03±0 


001 


11 


91±0 


001 


11 


68±0 


006 


11 


62±0 


005 


11 


62±0 


014 


11 


55±0 


038 





26 


-2 


71 


683 


4.8 




21 


04 


25 


63±0 


47 


50 


13 


40 


91±0 


04 


13 


33±0 


001 


12 


66±0 


001 


12 


42±0 


001 


12 


25±0 


006 


12 


22±0 


009 


12 


06±0 


039 


12 


27±0 


065 





17 


-2 


83 


684 


8.1 





''Magnitudes are de-reddened when extinction value A^ra quoted. 

determined from method outUned in §5. 
'^Slope a determined from IRAC photometry 
•^Running number from lHerbig &: Dahml l|2006h . 
°EW[Ho] from lHerbig fc Dahml l|2006h . 



Table 6. Assumed Reddening Law from 



Flaherty et al 



(120071). 



Waveband 


Ax/Aks 


J 


2.5 


H 


1.55 


K, 


1.0 


Band 1 (3.6 /im) 


0.64 


Band 2 (4.5 /^m) 


0.52 


Band 3 (5.8 /im) 


0.49 


Band 4 (8.0 urn) 


0.47 
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21 0.4'^ 30^ 03^ 30^ :p2^ 3p^ 

Right Ascension 
Center: R.A, 21 03 29,05 Dec +50 1.5 0.9.9 



Fig. 1. — Three-color 3.6, 5.8, and 24 /zm image of the young asymmetric cluster L988e, 
including regions f and c. The HAeBe stars LkHa 324 and 324SE are found in the L988e 
region. 
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Ri.qht Ascension 



Fig. 2. — Palomar Observatory Sky Survey II V-band image with the rectangular fields 
of view for the data presented in this work superimposed. Additionally, the embedded 
(blue) and exposed (red) regions of the cluster and the off-cluster fields (white) used when 
determining disk fractions are shown in this figure as circles. 
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21^ 0,4"^ 30^ DO^ 03^ 3Q^ 00^ 

Right Ascension (J2000) 



Fig. 3. — Three-color J, H, K ,,, image of the L9 88e region, with contours of ^^CO emission 
(50" resolution) overplotted ( feidge et all 120031 )). The ^^CO emission has a maximum inte- 
grated intensity of 13.75 K km s~^ and contours range from 10% to 90% of the maximum in 
increments of 20%. Squares designate objects classified via the {J)H,Ks,[3.6],[4:.5] methods 
or IRAC+MIPS methods and triangles designate objects classified via IRAC [3.6] and [4.5] 
and MIPS [24] photometry. Red designates identified protostars, green designates identified 
stars with disks, and blue designates identified stars with transition disks. Note the paucity 
of stars detected toward the molecular cloud in the image. 
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Fig. 4. — Linear grey scale SCUBA 850 fim emission (15" resolution), with contours of ^^CO 
(50" resolution) overplotted, with the same contours as in Fig. 3. Red designates identified 
protostars, green designates identified stars with disks, and blue designates identified stars 
with transition disks. 
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00 

Right Ascension 



03^ 30^ 



00' 

J.O Mags 



Colortior;; K-bcnd Ej<t|jictian 



Fig. 5. — Axs extinction map from WIRC data in grey scale, with the same ^^CO contours 
as in Fig. 3, and identified YSOs overplotted. The grey color scale extends from Aks = 
0.2 (white) to 3 (black) mag. In the region of the maximum CO intensity, SCUBA derived 
values are substituted for A^s- Identified YSOs are plotted with red squares denoting class 
Is, green squares denoting class lis, blue squares denoting transition disk objects, and blue 
x's denoting class Ills with Ha enaission . The YSO identified by filled symbols are Ha 
emission-line stars ( Herbig fc Dahm 20061 ). 
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1 20 ^ield/Class III ' 
100 



3 80 h ' 

o 60 h 
6 

z 40 

20 




Class II 



I'l 




5 



-2-10 1 
Dereddened 2-8yCxm SED slope 



2 



Fig. 6. — Histogram of the num b er of objects vs. de-reddened 2-8 /im slope for six young 
clusters included in iGutermuthI (120051 ). Class Ill/Field stars are well separated from class 
lis, and class Is are defined by consideration of distribution shapes. 
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Fig. 7. — Dereddened (i^s — [3.6])o vs. ([3.6] — [4. 5])o color-color diagram. Black dots represent 
either class III or field stars. A dashed blue line shows the separation between class I and 
II. Solid blue lines show the limits of regions clearly dominated by class I and II objects, 
as given by de-reddened 2-8 /im SED classifications, while the dashed green lines show the 
less restrictive region actually classified here, allowing for YSO variability. Triangles denote 
objects originally classified as class I based on dereddened {Kg — [3.6] )o vs. ([3.6] — [4.5] )o 
colors, but were reclassified as class II due to a ([5.8] — [24]) color of less than 4. 
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Fig. 8. — Histogram of the number of stars vs. de-reddened 3.6 /xm magnitudes for all 
objects classified as having an infrared excess from {J)H ,Ks,[3.6], [4.5] photometry. Note 
the separation between YSOs ([3.6] < 13.5) and extragalactic contaminates ([3.6] > 13.5). 
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Fig. 9. — NIR Color-Magnitude diagram for cluster member s (extragalact i c cont aminants 
removed) of the L988e region with the 1 Myr isochrone of iBaraffe et al.l (119981 ) overlaid 
in black. Solid lines are the reddening vectors for a 1 Mq object (top), 0.08 Mq object, 
the Hydrogen Burning Limit (middle), and 0.025 Mq object (bottom). Gray x's represent 
objects classified as class Ill/photospheric and black squares represent objects classified as 
YSOs by a J, H, Kg, [3.6], [4.5] classification. 
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Fig. 10. — Three-color [3.6], [4.5], [8.0] image of the L988e star formation region. Overplotted 
on the image are stars with disks (green boxes), transition disks (blue boxes), and protostars 
(red boxes). 
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Fig. 11. — De-reddened SED slope (from 3.6 fim to 8 fim) vs. EW[Ha] in A. Here red 
squares represent class Is, green squares class lis, blue squares transition disks and blue x's 
class Ills determined from (J),H ,Ks,[3.Q], [4.5] photometry. A vertical dashed line at 10 A 
illustrates the accepted single- valued equivalent width division between WTTS (< 10 A) 
and CTTS (> 10 A). Circled objects a re th ose that would yield a different CTTS/WTTS 
classification using the IWhite &: Basril (120031 ) criterion for CTTS/WTTS classification as a 
function of spectral type. Green circles imply CTTS rather than WTTS as determined from 
the single- valued equivalent width division, and blue circles imply WTTS rather than CTTS. 
Two horizontal dot-dashed lines at slopes of -1.8 and -2.56 delimit the 'anemic' disk region. 
Sources with thick disks fall a bove the -1.8 line and sources without disks fall below the 
-2.56 line. As demonstrated by iLada et al.l (120061 ) . sources with disks (and envelopes) have 
stronger Ha equivalent widths. 
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Fig. 12. — Grey scale SCUBA 850 /xm emission in greyscale, overplotted by contours of the 
surface density distributions of classified cluster members. The locations of LkHa 324 and 
324SE are given by an astrisk and an x, respectively. The most compact distribution, the 
embedded population, with a relatively high percentage of class Is, is concentrated on the 
dust maximum. To the east of the dust ridge, there is a second, les s compact conceritratio n, 
of primarily class lis, associated with the emission-line population (IHerbig fc Dahnj|2006l ). 
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Fig. 13. — Top Left: histograms of i^'g-band detections in 1.5' diameter circle centered 
on an off-cloud region (solid line) and on the exposed cluster (dashed line). Top Right: 
histograms of remaining 22 objects after subtracting the off-cloud sample from the exposed 
cluster sample (black, solid line) and the known population with disks (red, dot-dashed 
line). Bottom Left: histograms of i^^-band detections in 1.5' diameter circle centered on 
an off-cloud region (sohd line) and on the embedded cluster (dashed line). Bottom Right: 
histograms of remaining 19 objects after subtracting the off-cloud sample from the exposed 
cluster sample (black, solid line) and the known population with disks (red, dot-dashed line). 



